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Abstract 

We study the renormalisation of SU (W) gauge theories on general anisotropic lattices, 
to one-loop order in perturbation theory, employing the background field method. The 
results are then applied in the context of two different approaches to hadronic high-energy 
scattering. In the context of the Euclidean nonperturbative approach to soft high-energy 
scattering based on Wilson loops, we refine the nonperturbative justification of the analytic 
continuation relations of the relevant Wilson-loop correlators, required to obtain physical 
results. In the context of longitudinally rescaled actions, we study the consequences of one- 
loop corrections on the relation between the SU(Nc) gauge theory and its effective description 
in terms of two-dimensional principal chiral models. 


1 Introduction 

Anisotropic lattices are a standard tool in modern lattice calculations, and have been used in the 
study of a large variety of problems, ranging from glueball [T] and light-hadron [2] spectroscopy 
to properties of QCD at finite temperature [311!. Numerical calculations in four dimensions 
usually employ lattices with 3-1-1 anisotropy, i.e., only one of the lattice spacings is different 
from the others, while more general anisotropy classes have received much less attention [5], due 
to the increasing difficulty in the scale setting procedure. Indeed, for anisotropy classes other 
than 3-|-l, one needs to appropriately tune the action in order to recover Lorentz invariance in 
the continuum, already at the pure-gauge theory level. A better understanding of these more 
general anisotropy classes would be useful, since they provide a more flexible setting for varying 
length scales independently in different directions. This would allow, for example, to enlarge 
the range of momenta accessible to lattice calculation at a reasonable computational cost, by 
improving the resolution only in a single spatial direction [^. 

Anisotropic lattices provide, quite obviously, the natural setting for the nonperturbative 
study of anisotropic systems, also beyond numerical applications. An interesting case is that of 
longitudinally rescaled actions, which in recent years have been considered in the context of high- 
energy scattering in QCD [3 [Z1131313 HD IE] . The basic idea of Refs. [3 0 1313 is to perform a 
rescaling of the longitudinal directions, which appear highly Lorentz-contracted in a high-energy 
scattering process, in order to derive an effective action starting from QCD. In Refs. I30S 
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[9] only the classically rescaled action was considered, while the important effect of quantum 
corrections was studied later in Refs. [ini ttu IE], in the framework of Wilsonian anisotropic 
renormalisation in the continuum. In this context, the use of a gauge-invariant, anisotropic 
lattice regularisation could lead to more insight in the structure of quantum corrections. The 
relevant anisotropy class here is 2+2, with different lattice spacings in the longitudinal and in 
the transverse plane. This is also the case considered in Ref. [5], although for different purposes. 

In Ref. [13] a classical anisotropic rescaling of the functional integral has been used to justify, 
on nonperturbative grounds, the analytic continuation from Euclidean to Minkowski space, 
required to obtain physical results in the Euclidean formulation [HI ESI (III dzi Eg [H] of the 
nonperturbative approach to soft high-energy scattering [2ni|2ll[22l|23l|231[25l|26l|27j. This 
approach has been recently used in Ref. [28] to obtain a theoretical estimate of the leading energy 
dependence of hadronic total cross sections, resulting in fair agreement with experiments. As the 
analytic continuation plays a key role in this approach, it is important to establish its correctness 
going beyond the formal argument of Ref. m. which, as we have said above, is based only on a 
classical rescaling of the QCD action. To this end, quantum corrections to the effective action 
must be included to prove that the necessary analyticity requirements are actually fulfilled. The 
relevant anisotropy class in this case is 2-|-l-|-l, with different lattice spacings in the transverse 
plane and in the two longitudinal directions. 

The purpose of this paper is to perform the renormalisation of a SU{Nc) gauge theory 
regularised on a general anisotropic lattice, and to apply the results in the study of hadronic high- 
energy scattering through the approaches mentioned above. To avoid the complications related 
to the introduction of fermions on the lattice, we work here in the quenched approximation, i.e., 
pure gauge theory. 

The plan of the paper is the following. In Section [2] we study renormalisation for a general 
anisotropic lattice regularisation, using the background field method on the lattice [29l [3011311 [32l 
IMIIMIIMIIMIIST] . In Section[3]we use the results in the context of the nonperturbative approach 
to soft high-energy scattering of Refs. [20l|2ll[22l|23l[Sl|25l[26ll27], refining the argument of 
Ref. [13] on the possibility of performing analytic continuation to Euclidean space. In Section 
H we discuss the longitudinally rescaled actions of Refs. P EllSl E] EE] EH ES], focussing on the 
representation of the SU{Nc) gauge theory as a set of coupled two-dimensional principal chiral 
models. Finally, Section [5] contains our conclusions and prospects. Some technical details are 
discussed in the Appendices. 


2 Anisotropic renormalisation 

Our aim is to renormalise the Euclidean SU{Nc) gauge theory regularised on a 4D orthogonal 
anisotropic lattice. More precisely, lattice points are located at x = A; where fi are four 

orthogonal unit vectors, and the physical coordinates = x^(n) in Euclidean space are x^(n) = 
a^n^, Ufj, G Z. Here = a/X^ is the lattice spacing in direction fj., with the dimensionless 
anisotropy parameters G being the inverse ratios of to a reference length scale o0 

^The five parameters a and {A,i} are obviously redundant, and some condition has to be imposed on {A^} to 
remove this redundancy. This notation is however convenient, as we treat all the directions on the same footing. 
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Consider the following Wilson-like action, 


1 


-5iaT = /3 E E 1 E E C>.uV,.{n), 


fj,iy \' '"' J j f-' '-' fJLU ' fJjiy \ ‘ “ J 7 ( 2 . 1 ) 

n ^<u \ c ' n 

where U^u{n) are the usual plaquette variables built up with the link variables U^{n) E SU{Nc), 

U^.{n) = U^{n)U,{n + A)C/^(n + v)Ul{n ), (2.2) 

j3 = 2Ncjg‘^ with g the coupling constant, and g i', are the plaquette coefficients^ 

4 


J, J-^ = \{X^. 


(2.3) 


a=l 


It is straightforward to show that Eq. (12.11) yields the correct naive continuum limit upon iden¬ 
tification of the continuum, physical gauge fields A^{x) through 

U^{n) = ^ (2.4) 


as appropriate for an anisotropic lattice. The choice of plaquette coefficients C^u is easily 
understood by noticing that J is just the volume of an elementary cell, so that J is the 
Jacobian for the change of variables from isotropic to anisotropic coordinates, while = 

{X^Xy)~^a? is the area of the faces of an elementary cell lying in the gv plane. 

As is well known, divergencies appear in the continuum limit when taking into account quan¬ 
tum corrections. These divergencies need to be subtracted through a suitable renormalisation of 
the couplings in order to obtain a finite continuum theory. On the isotropic lattice, the symme¬ 
try under the unbroken hypercubic subgroup of 0(4) guarantees that all the plaquette terms in 
the action need to be renormalised in the same way, so that a single redefinition of g is sufficient 
to reabsorb the divergencies. The form of the action is therefore unchanged, and one recovers 
the full 0(4) invariance in the continuum limit. 

On a general anisotropic lattice this residual symmetry is broken, except for reflections 
through lattice hyperplanes, and so in general different terms will require a different renormal¬ 
isation. Since there are six different plaquette terms and only four lattice spacings, it will not 
be possible in the general case to reabsorb completely the quantum corrections into a redefini¬ 
tion of A^, keeping at the same time the same form of the tree-level action [5]. In turn, this 
implies that the continuum limit of Eq. ()2.1I1 cannot be made into an 0(4)-invariant theory by 
an appropriate, simple rescaling of the lattice spacings, since in the general case one will still 
find different coefficients for the six continuum field-strength terms. To recover 0(4) invariance 
one must ensure that these coefficients are equal, and this requires that we take the action to 
be of the more general form 

5’iat = E E = E E 1 ^> (2-5) 

n fj.<u \ c ' n ^<u 

^For definiteness, we define also = 0. 
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where the couplings (3^i, = = PfiuW have to be properly tuned to yield a finite, 0(4)- 

invariant theory in the continuum limit. 

The need for tuning comes, as we have said, from the fact that there are in general more 
couplings than anisotropy parameters. It is however easy to show that one has to tune at most 
only two combinations of the couplings to achieve restoration of 0(4) invariance in the continuum 
limit, while the other four independent combinations can be interpreted as the coupling fixing 
the overall lattice scale, and renormalisations of the anisotropies a^ja^ = To see this, 

let us remove the redundancy in the set {A^} by imposing the symmetric condition ]))[^ = 1, 

thus defining a in terms of the volume of an elementary cell. Any other equivalent choice (i.e., 
giving the same a^) is obtained by a simple global rescaling of {A^} and of a. The six plaquette 
terms can be grouped in pairs of “complementary” fj,iy and fLu plaquettes, i.e., {Ui 2 ,U 3 i), etc., 
which we denote as (/ij/|/uz/) = (12|34), (13|24), (14|23). It is also easily noticed that Cf^^, = 
so that Cpp = C~^. This suggests to parameterise /3^i, as follows, 

Pfiiy = /3 ~~ ■ ( 2 - 6 ) 


As there are two redundant parameters, we choose to fix = 1, so that our condition 

on {A^} is not renormalised H and n(/ii/|^p) are 

unambiguously defined and can be obtained from as follows, 






Z, 








(2.7) 


This makes it explicit that the restoration or not of 0(4) invariance in the continuum depends 
only on the values of the ratios of the couplings (3^i,. Defining now the bare anisotropy parameters 
A^ = and the bare plaquette coefficients Cj^^, = C'^j/(A'®), one can rewrite Eq. (12.5p as 


5iat = /3 ^ Z(^,|^,)[0^,P^,(n) + 0|,P^,(n)] . (2.8) 


This equation shows that to obtain an 0(4)-invariant theory in the continuum limit, one can 
choose freely A^ (up to a constraint to remove the redundancy), and then tune only the two 
independent ratios of to the appropriate values. The physical anisotropy parameters A^ 

are related to the bare ones through the renormalisation A^ = z~^Xj^, and can be measured ex 
post. 

Using the parameterisation Eq. (1231), it is possible to set up a rather simple nonperturbative 
scheme to achieve restoration of 0(4) invariance in the continuum, for an arbitrary choice of bare 
anisotropy parameters. The basic idea is to impose that the string tension, determined from the 
asymptotic behaviour of large rectangular T x R on-axis Wilson loops Wa/s ~ exp{— 
is the same for all pairs of directions a, (3. Denoting with a the physical (dimensionful) string 

®Any other choice is of course allowed. If, for example, the scale a is defined to be one of the lattice spacings 
by choosing = 1 for some fi, then it is convenient to choose = 1. The new values of Zi, are obtained from 
those corresponding to the symmetric condition by replacing Zi^ —>■ while and /3 are unaffected. 


4 








tension, this amounts to impose Xa^/3^a/3 = o?a, for all pairs of different a, j5. Multiplying the 
relations for dap and its “complementary” dap-, one obtains the following consistency conditions, 

dl2d34 = (Ti3d24 = di4(T23 , (2.9) 

which have to be imposed to recover 0(4) invariance. This can be done without any prior 
knowledge of the physical Aq, and requires only to properly tune two of the coefficients 
(the third one being constrained by our choice n(/^iy|gi 7 ) ^{iiu\p.p) — !)• Having done this, the 
anisotropies can then be obtained from the ratio X^/Xu = for any a / Imposing 

](([^ = 1 one can explicitly determine all A^’s, and set the lattice scale a from the relation 

a^cj^ = (Ti 2 d' 34 . While the string tension is known not to be the best observable for setting the 
physical scale, nevertheless it could be useful for the tuning, as it can be determined to high 
precision by means of multilevel algorithms [38]. It is worth mentioning that the tuning of two 
parameters is only required when all the lattice spacings are different: if at least a pair of lattice 
spacings are equal, one easily sees that only one parameter has to be tuned0 

2.1 Background field method 

From the discussion above, we see that our task is to find the relations among the couplings 
that will lead to an 0(4)-invariant theory in the continuum limit. We will study this problem to 
lowest order in perturbation theory, making use of the background field method [29113(11 132] 

on the lattice |33l [Ml Ea ESI ETj. The advantage of this method is that it allows to keep an 
exact gauge invariance on the lattice after gauge fixing, which greatly simplifies the calculations. 
A full account on the background field method can be found elsewhere [391BO] . Here we briefly 
recall the main points of the method to fix the notation. 

(c) 

The first step is to introduce a background field in the gauge action as follows, 

Sbf[U^'^\v] = , ( 2 . 10 ) 

where we now denote with V the gauge links, to be integrated over with the usual Haar measure. 
As a consequence of the gauge invariance of 5iat, the action 5 bf is invariant under the background 
gauge transformation 

t/W«(n) = G{n)Ujf\n)G\n + fi), V^{n) = G{n)V^{n)G\n), (2.11) 

with G{n) G SU{Nc), as well as under the following gauge transformation of V alone, 

H^(n) ^ G{n)V^in)Ulf\n)G\n + . (2.12) 

The integration measure is also obviously invariant under the transformations Eqs. (|2.1ip and 
(j2.12p . One then proceeds to set up perturbation theory in the usual way, setting 

Vf,{n) = , Hjf)(n) = ^ (2.13) 

^ In the 3+1 case, where a single lattice spacing differs from the others, there are only two kinds of plaquette 
terms and so only two independent lattice string tensions. In this case there is thus no consistency condition to 
be satisfied and no tuning is needed, as is well known. 
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wher^ Quin) = q'^{n)t°‘ and B^{n) = with t“ the generators of SU{Nc) in the funda¬ 
mental representation, a = 1,... ,N^ — 1, and tr{t“t^} = One then changes variables of 

integration to q, expressing the Jacobian as a contribution 5meas[<?] to the action. Notice that 
powers of g and a are chosen so that q is dimensionless, while B has dimensions of mass. This 
distinction is convenient for book-keeping purposes [55] . 

Under the transformation Eq. ()2.1ip . the background field B transforms as a gauge field, 
while the “quantum” field q transforms as a matter field in the adjoint representation. The 
symmetry under the gauge transformation Eq. (12.121) requires to impose a gauge condition on q 
to define the corresponding propagator. This is done d la Faddeev-Popov, adding a gauge-fixing 
term to the action, together with the corresponding ghost term. The key point is that there 
is an appropriate choice of gauge, called the background field gauge, for which the gauge-fixing 
and the ghost terms are invariant under the background gauge transformation Eq. (12.lip . This 
gauge-fixing term is [29l[3Ql[37] 


Sg.i.[B,q\= ’ 


where are the lattice background covariant differences, 


Dfi fin) = 


Ulfi\n)f{n -h jl)Ulfi'^\n) - /(n 

- fi)fin - fi)Ulfi\n - fi) - /(n) 


(2.14) 


(2.15) 


in which a factor is also included for convenience. The usual lattice differences are 

(c) 

obtained setting = 1 in the expressions above, where 1 denotes the unit matrix. The 
corresponding ghost term is 


5’ghost[-B,g,c,c] = 2J^^tr |[T)+c(n)] M ^ (^■^q^in)'j D+ + iAd (^^q^{n)'j c(n)} , (2.16) 




where c = c = with c“, c“ independent Grassmann variables, and where 

1 _ e-iAd(X) 

Mix) = , Ad(X)y ^ [A, U]. (2.17) 

It is straightforward to prove invariance of these two terms under the background gauge trans¬ 
formation, Eq. (|2.11l) . supplemented by the transformation laws for the ghost fields, 

c^(n) = G(n)c(n)G^(n), c‘^(n) = G(n)c(n)G^(n). 

Expanding Eq. (|2.16l) up to Oig^), one finds 

Sghost[B,q,c,c] = 2J ^tr {[T>+c(n)][T>+c(n)]} +Oig) 

= 2j7^tr {c(n) DfiD^cin)} +Oig) = S’°host[-S, c, c] -b 0 ( 5 ), 

®Here and in the following, the sum over repeated colour indices is understood. 


(2.18) 


(2.19) 
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where we have used “integration by parts” on a lattice (infinite or with periodic boundary 
conditions), 

^ tr {/(n)[L>"c/(n)]} . (2.20) 

n n 

The starting point for the perturbative analysis is the generating functional 

Z[B, J, fi,v] = J 

Stot[B,q,c,c] = SBF[B,q] + 5*meas [ 9 ] T 'S'g.f. [.B, qj + S’ghostfB, q, c, c], 

where with a small abuse of notation we have written 5BF[B,q'] = and we have 

added source terms for the various fields. Here J = J{n) = and J-q = J2n ti 

and similarly for the other terms. A Legendre transform gives the effective action (generating 
functional for IPI graphs), 


r[B,Q,C,C] = -W[B,J,f),r,] + J-Q + f)-C + C-g, 


( 2 . 22 ) 


where the classical fields Q, C and C are defined as 


Qlin) = 


5J“(n) 


C“(n) = 


9fL[B, J, r/, 77 ] 
(n) 


C“(n) = 


dW[B,J,f),r]\ 

dr]^{n) 


(2.23) 


i.e., they are the expectation values of the quantum fields for prescribed values of B and of the 
sources. 

Defining a background gauge transformation for the classical fields, imposing that they trans¬ 
form as the corresponding quantum fields, Eqs. (12.lip and (j2.18p . leads finally to the identity 

r[B^, C^, C^] = r[B, Q, (7, C] (2.24) 


for the effective action. This is the key relation that allows us to simplify the calculations. 
Indeed, setting 5'efr[B] = r[B, 0,0,0] — r[0, 0,0,0], as a consequence of the bacl^round gauge 
invariance, of the discrete symmetries of the action (translations and reflection^, and of the 
locality of divergencies, to one-loop accuracy and to lowest order in perturbation theory we are 
guaranteed to find in the continuum limit 


lim 

a^O 


= 5 E / 


d^x 


A 


_EL -K 
2Nr 




-|- (non-local finite terms) -|- 0{g ^), 


(2.26) 


where = Ky^ = K^y{a,X) is 0{g^), and where B^y = d^By - dyB^ -(- i[B^,By\ is the field 
strength for the continuum background field B^{x), B^(x(n)) = B{n). For our purposes it is 

^Reflections lie, act as follows on the coordinates, ITan^ = for fj, ^ a, Tlaria = —ria- The corresponding 
transformation laws for B and q are the following, 


-Bc,{Ucn-&), ^ ^ \ -Ui‘''^\n^n-a)q^{n^n-a)Ut\Uc.n-a). 


( 2 . 25 ) 
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therefore sufficient to compute the two-point function of the background field to have enough 
information to renormalise the theory and impose 0(4) invariance. To one-loop accuracy it is 
enough to set 


_|_ 

2Nc “ g2 + 2Nc 


— K — _ 

i/r 


(2.27) 


where gr is the renormalised, A-independent coupling. 

We notice that Eq. (12.51) , with the couplings chosen according to Eq. (12.271) , can be interpreted 
in two ways. Under the identification t/^(n) = with Xa = CLana, it leads in the 

continuum to the renormalised, isotropic action for the gauge fields A^{x), for which it provides 
an appropriate lattice discretisation. On the other hand, identifying U^{n) = with 

Ua = aria, in the continuum limit one obtains the following renormalised anisotropic action, 




1 




(2.28) 


with <l>^jy the usual field-strength tensor for (j)^, for which Eq. ()2.5p provides therefore a lattice 
discretisation. This is the form of the action obtained by classically rescaling coordinates and 
fields in the Yang-Mills action, discussed in Refs. laiTlElIlKIS]. 


2.2 One-loop calculation 

To compute it is enough to expand the action to order 0{g^), which in turn means expanding 
the gauge action up to second order in q. Contributions from S'meas are at least 0{g‘^) and can 
be ignored. Let us expand the action S'bf + 5'g.f. in powers of g, 

SriY[B,q] + 5g.f.[i?,g] = Sc[B] + S^i[B,q] + S^ 2 [B,q] + ... , (2.29) 

where S'c[i?] = S'bf[.B, 0] is the classical action, S'gi[R,g] is linear in q, Sg 2 [B,q] is quadratic and 
so on, and set 

S 2 [B, q, c, c] = Sg 2 [B, q] + c, c] 

n^m 

A straightforward calculation then shows that 

Ses[B]\.^i 0 ^ = Sbf[B, 0] + - log - log . (2.31) 

etn Jlo(gO) uri , j 2 ^ detn[0] detn[0] 

Terms linear in q play no role and can be ignored0 Eq. (j2.3ip can be conveniently written as 

(2.32) 

^These terms are usually discarded by requiring B to satisfy the equations of motion, but this is actually not 
necessary. 











where c, c] = S 2 [ 0 ,q,c,c\ is the free action with no background field, and (.. .)o denotes 

the corresponding expectation value, 


{0[B,q,c,c])o = / VqVcVce 




0[B,q,c, c] , 


Zfree = / PgPcPc 


(2.33) 


For future utility, we define the connected correlation function (0i02)oc = {Oi 02 )o — {Oi){ 02 ) 0 - 
Since we are interested only in the two-point function for B, only terms up to 0{B^) will be 
kept in 82 - 


2.2.1 The quadratic action 

The gauge action in a background field can be conveniently written as follows, 

5BF[S,g] = (l - ±Retr {V(n)C/W(n)}] , 

n,ij,<u \ c / 

where the “quantum” and the “background” plaquette are given respectively by 

C/W(n) ^ . 


(2.34) 


(2.35) 


A standard application of the Baker-Campbell-Hausdorff formula gives 

C/W(n) = exp |i^/^,(n) + 0{a^BdB, a\dBf) + 0{a^B^)^ , (2.36) 

witlil 

= a-i {A+B, - A+B^) + i[B^, B,], (2.37) 

which in the continuum limit reduces to the usual field strength tensor for the background field 
For we have instead 


where 


V^uin) = exp < ig 




[Fiivin] + gR^,y{n)\ +0{g'^) 


Ffj.1/ — D^ qv Du 9ai) 

= l^T^[D~llqu,D^q^]+i[q^,qu\, 

D?u = \ ’ 


(2.38) 


(2.39) 


®In the following equations we will sometimes drop the dependence on the lattice site n to make the expressions 
more readable. 

® In principle, also the higher-order terms of order 0{a^BdB, a^{dE)^) appearing in Eq. (12.361) could contribute 
to the two-point function in the continuum. This however is not the case, as we will see below (see footnotes [iQ] 
and fTTI) . 
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(2.40) 


and = R^i^} + i?®. Expanding up to quadratic terms in B and q we find 
5BF[-B,g] = Sc[B] + Sq[B,q] + (linear in q) + 0{q^), 
where Sc is the classical action, already defined above, 


Sc = Ja^ -tr/^ (n) 

2Nc 2 ^ a- 




while the “quantum” piece Sg is given by 



Sq[B,q] =JY^ ^tT{F^^{n) + 2a‘^Rf,y{n)f^y{n)} 


1 a‘^ 

4 (A^A,)2 


tr{i^;i(n)4^^(n)} . 


(2.41) 


(2.42) 


The gauge-fixing term is quadratic in q, and can be conveniently rearranged as follows, 

‘^g.f. = *5^'f^. + <52. + *5T', (2.43) 

where 

^ X] {DyqtM{n)D+qy{n)} , = Jti: {R^f}Kn)ftMu{n)^ , 

RI}u = i (^[Qt^,Qu] + Y[q,i,D^qu] - Y[qu,D^qf,] + j^[D+qf„D+qy] 

St' = ja^ Yl tri?!S(n) 

^iiu ~ r)\ \ ~^^^qv + qu\[f^lVj 9m 4" 9m] • 

Finally, the ghost term is independent of q to 0{g^). Putting all the terms together, one obtains 
for the quadratic lattice action 

52 = 5e + + Sk +Sb +St + St' , (2.45) 


(2.44) 


where the terms have been grouped so that each quantity in the equation above is separately 
invariant under a background gauge transformation [36] . Here = •S'gj®® ^ , with 

^gfuon = ^ E [A+9.(n)]' , 5fr®®g, = 2 J y; tr {c(n) A;A+c(n)} , (2.46) 


being the free actions for gluons and ghosts, respectively, in terms of which the propagators are 
defined, while the interaction terms are given by 


'^g’fuon = Sq + S^^l - S’^^®®n = ^ E tr {[-D+g,.(n)]2 - [A+qy{n)f} , 
S^host = S-g^ost - 5’^host = 2jr E tr {c(n) [D^D; - A;A+] c(n)} . 


(2.47) 
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Moreover, extra vertices come from the terms 


Sk = {R^i^}{n)f^u{n)] + sf} = Ja? tr {[/^^(n) + -R{,l)(n)]/^^(n)} , 

5b = Ja^ Yl (4-(n)/M.(n)} , = -Ja^ Y m \ {^/i(n)/A(n)} ■ 

Explicitly, we have for 5 a and Sb the expression^ 


Sa = Ja^Y {^/,t.(n)/^,,(n)} , 

— 2i ^2, Qi^] + [qfj_ , qiy] — -j— [qiy , q^] 

Sb = Ja^ Y ’ 


2 A^A/ ^ 


[DU,,DtqY, 


B^u — * I \ 
A 


(2.48) 


(2.49) 


Notice that the terms 5 a and 5 t' are odd in a given component of the gluon field, while 
the other terms are even. Since the propagator is diagonal, this implies [33l ESI ES] tha10 
(*S'a)o = (*S't')o = 0, and also that (5"|*);^ jj 5 a)oc = (*S'a-S'b)oc = 0. 


2.2.2 The effective action 

Expanding now Eq. (I2.32p up to terms quadratic in the background field, we obtain the following 
expression for 5eff[i?], 


Reii\o{B^),O{g0) 



- 2(('^b)^)oc + (*S't)o + ((5'b)o - (5'gfuon5'B)oc) 


= 5e + A5g + AS A + A5b + A5t + A5gB • 


(2.50) 


Here we have taken into account the remarks after Eq. (12.491) , and the fact that in four dimensions 
the ghost contribution exactly cancels half of the gluon contribution from [331 ESI ESI ES] • 

Terms have been grouped so that each contribution is separately gauge-invariant [36] . 

The evaluation of the various terms is performed generalising the techniques developed in [3S| 
to the anisotropic case. Since such a generalisation is straightforward, here we simply list the 

In these quantities one should in principle include also the higher-order terms mentioned above in footnote 
[9l by properly redefining 

This clearly remains true also if higher-order terms neglected in Eq. (12. ,361) are included in the definition 
of Z/ii/, see footnotes [9] and [To] Since S'gYuon is 0{B), the only contribution of a higher-order term which should 
still be considered is that of the 0{a^BdB) term in Eq. (12.361) to (S'b)o oc {{Bf^Ao^}f^v\ however, (global) 
background gauge invariance implies that tr {(i3^i/)ot“} oc 5“, and so higher-order terms can be safely ignored. 


11 













term 


AK, 


fll' 


AS, 

ASa 

ASb 


AS^ 


Nr 


3(47r) 




- W - - 5;.W + 5W1 


4A^c 

(47r)' 


-7 + log 


{aAiy 


+ 2 


+ Nr - 2g^{x) - 2g,ix) + Ag{x)] 


4 




AJ A57 A2 AJ 


- 1 

ASt " 


2 1 r^^(A) ^ z,(A) 


2A, 


A2 A2 , 


Table 1: Contribution of the various terms in Eq. (j2.50p to Kfj_iy in the effective action, Eq. (j2.26p . 


results, giving in Tab.[I]the contribution AK^y of each term to the quantity K^y [see Eqs. (|2.26p 
and (12.271) ] in front of 1/2/ d^xiiF^y{x). The relevant technical details can be found in the 
appendix of Ref. [36]. Summing up, one obtains 

K^y{a, X) = + lC^y{X) = /3o log + /C^,(A), (2.51) 


with M a mass scale which sets the renormalisation point, /Iq the first coefficient of the Yang- 
Mills /3-function im 021 113], 


/3o 


11 Nr 
3 (47r)2 ’ 


(2.52) 


and with KL^y finite, a-independent coefficients. 


a:^.(a)= 


3 (47r)2 

N, ^ 


r 641 


1 

+ 

CO 1 
CO 1 

+ Nr 


^g,yix)-^{g,ix) + gy){x) + ^g{x) 


+ 


Z(A) 


1 1 


^/a(A) ^.(A) 


A^ 




+ 


N^-l 

2Nr 


Z^iX) ^ Zy{X) 


XI 


XI , 


(2.53) 


where 7 ~ 0.5772 is the Euler-Mascheroni constant, and 0^,^, 0^, g, Z^ and Z are functions 
of {A^} defined in terms of integrals involving the modified Bessel functions of the first kind. 
Their precise form is not needed for the analysis of the present Section, and can be found in 
Appendix lAl Eqs. (lA.ip and (IA.3E 

To renormalise the theory and recover 0(4) invariance in the continuum limit it is enough 
to set 

1 . . 1 . 1 <3/3;.. 


1 


r 


gKM) 


+ /3o log 


{aAiy 


= Po log 


{aAy 


2Nr 


= JC 




(2.54) 
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Here A = M exp{—1/(2/305^ (M))} is a renormalisation-group-invariant mass scale, whose value 
can be determined by comparing lattice results with experiments. Since a shift <5/3^;^ —>■ + 13 

can always be reabsorbed in a redefinition of g, any choice satisfying the set of conditions <5/3^;^ — 
^(3pa = ICpi, — Kpa will actually lead to restoration of 0(4) invariance at one-loop accuracy!^ 
As we show in Appendix!^ under a global rescaling ^ Z and Z^ get a factor 
and Qp are unchanged, and Q ^ Q + ^^logC^, so that overall XZpy —)• l+py + /3ologC^- Since 
the additive term can be cancelled by a —>■ this means that the couplings /3^jy depend on 
a and Xp only through the combinations provided by the lattice spacings a^, as they should. 
As we have already remarked, to avoid redundancy one has to impose a condition on the A^’s, 
like, e.g., setting Aq = 1 for some a, so using one of the lattice spacings as reference length, or 
imposing the symmetric condition Aq, = 1, thus using the volume of the elementary cell to 
define a. 

We have compared our results with the ones available in the literature for the isotropic |33l 
[Ml [36], 3-1-1 [Ml 05 SSI [5] and 2-|-2 |5| anisotropic cases F5 In particular, we have successfully 
checked that in the isotropic case we recover the result of [Mj, and we have compared the 
differences of 5/3pu with the ones reported in Ref. [5] for the 3-1-1 and 2-|-2 cases. While there 
is full agreement for the 3-1-1 case, we found a discrepancy in the analytic expression of one of 
the two independent differences in the 2-|-2 case0 On the other hand, the numerical values 
also reported in Ref. |S| agree with ours. It has to be noted that the analytic result reported 
in Ref. |5] for that difference does not vanish when there is no anisotropy, as it should, so most 
likely it contains some misprint. 

For future utility, we report the lowest-order approximation for the expectation value {Vpi,) 

i-S-q 

of the plaquette terms. Setting Up{n) = e and expanding in g, one finds 

{Vp.) = ^X^(tr F2 ,)o + 0{g ^), (2.55) 

where Fpy = A+gj, — A'^qp [see Eq. (j2.39|) ]. and (.. .)o has been defined in Eq. (12.331) . A 
straightforward calculation yields 


{Vpu) = g 


-1 
2Nr 


ZpW , Z,{X) 
XI 






+ 0{g^ 


(2.56) 


^^More generally, it is the ratios fipa that will be constrained by the request of restoration of 0(4) invariance, 
see the discussion in Section [^ 

^®In the 3-1-1 anisotropy class one lattice spacing differs from the other three, e.g., + \\ = = As, while in 

the 2-1-2 class the lattice spacings are equal pairwise, e.g., A 4 = Ai 7 ^ A2 = A3. 

In the notation of Ref. [5] , the discrepancy is in , in particular in the coefficients of the quantities 

B^{2, 1 ) and B^{2, 1 , 1 ), for which we find respectively and + ^). 
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3 Analytic continuation in the nonperturbative approach to soft 
high-energy scattering 

In this Section we use the results of Section [2] in the context of the nonperturbative approach 
to soft high-energy scattering. After a brief review of this approach (the interested reader 
can confer Refs. [20l [211 1221 (231 [Ml 123 [23 [2Z] for a more detailed discussion), we discuss its 
formulation on a Euclidean anisotropic lattice, and we refine the arguments of Ref. m on the 
analytic continuation back to Minkowski spacetime. 


3.1 Euclidean approach to soft high-energy scattering 

Soft high-energy scattering is characterised by small transferred momentum squared t, |t| < 
1 GeV^, and very large total center-of-mass energy squared s, s S> 1 GeV^. In the approach 
of Ref. [23; hadronic scattering amplitudes in the soft high-energy regime can be obtained 
from partonic scattering amplitudes after folding with appropriate hadronic wave functions. In 
particular, for meson-meson scattering the basic quantity is the scattering amplitude of two 
colourless transverse dipoles, which in the soft high-energy regime is given in impact-parameter 
space by the correlation function of two rectangular Minkowskian Wilson loops [211 [22]. These 
Wilson loops are computed on the paths described by the classical trajectories of the dipoles, 
so forming a large hyperbolic angle x longitudinal plane, and are cut at proper times 

±r for infrared regularisation purposes [3- In turn, their (Minkowskian) correlation function is 
obtained after analytic continuation in the angular variable and in the length of the loops from 
the correlation function of two Euclidean Wilson loops of length 2T forming an angle 9 in the 
longitudinal Euclidean plane na [a [13 [13 ca iiEiiini- This approach can be generalised to 
describe scattering processes involving baryons [laaaEaaaaia. As the construc¬ 
tions and the arguments of this Section are easily adapted to this case, we restrict the discussion 
to meson-meson (dipole-dipole) scattering for simplicity. 

The relevant Euclidean correlator is given bjcj 

gE{e, T; Tx; /i; /2) = - 1, (3.1) 


where (.. .)e denotes the expectation value in the sense of the Euclidean functional integral, 
z± is the impact-parameter distance between the dipoles, and Ri± and fi are the transverse 
size of the dipoles and the longitudinal momentum fraction of the quarks in the two mesons, 
respectively (“dipole variables”). The Wilson loops W) 2 computed on the following paths 
(see Fig. [I]), 

= ±^ 1 "^ + ^ + = ±^ 1 "^ + df ^ 2 ) 

CF^ : X^2i''') = -I- /^R2 = '^U2T + df , 

with r E [—T, r], and closed by straight-line paths in the transverse plane at r = TT. The 
four-vectors ui ^2 are chosen to be ui ^2 = (± sin 0 _l, cos |), 9 being the angle formed by the 
two trajectories, i.e., ui ■ U 2 = cos9. Moreover, Ri = (0,Rj_L,0), z = (0,Tj_,0) and f^ = l — fi, 

^®Here and in the following we denote by wx a two-dimensional vector in the transverse plane. 
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/j = —fi, with fi G [0,1]. The Minkowskian correlation function is obtained from Eq. (j3.ip by 
means of analytic continuation as follows [ISE], 

Gm{x,T; z±-, Ri±, fi] R2±, f2) = GE{-ix,'iT-, z±; Ri±, fi] R2±, f2) • (3.3) 

Physical amplitudes are finally obtained from Qm in the limit T —>■ oo, and for asymptotically 
large y logs. It is worth mentioning that combining Eq. (13.3p with the 0(4) symmetry of the 
Euclidean theory one obtains the following crossing-symmetry relations [mug], 

GMix,T-, Z±-, Ri±, fr, -R2±,1 - /2) = GiviiiT^ - X,T; Z±-, Ri±, fi; R 2 ±, f 2 ), (3.4) 

which allow us to relate the scattering amplitudes in the direct (meson-meson) and crossed 
(meson-antimeson) channels. 

The analytic continuation relation, Eq. (|3.3I) . has allowed studies of the correlators through 
nonperturbative Euclidean techniques [28l HU 081 09l HOI EH EZl EH EH EH ES] . For a brief 
review of the older results and a comparison to lattice data cf. [SHEIESI. 

3.2 Anisotropic lattice formalism 

It is well known that the functional integral needs to be regularised to become a well-defined 
mathematical object. Furthermore, the analytic continuation relation Eq. (j3.3p is meaningful 
only if a sufficiently wide analyticity domain exists. The first issue can be dealt with by discretis- 
ing the theory on a lattice, so that the relevant Wilson loop correlator can then be computed 
nonperturbatively, for example by means of numerical simulations, using off-axis operators to 
approximate the continuum Wilson loops. Numerical simulations using an isotropic lattice have 
been reported in Refs. [SHEIIES]. Unfortunately, only a discrete set of angles is accessible in 
this case; furthermore, for each angle one has to use a different off-axis Wilson loop, which makes 
the angular dependence even less analytically controllable. Since our purpose here is to study 
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the analytic dependence on 6 and T, it is more convenient to use an appropriate anisotropic 
lattice keeping fixed the Wilson-loop operator, which allows us to expose the dependence on 
the relevant variables in the action. In this way we make the functional integral a well-defined 
object, and at the same time we can study the analyticity domain of the correlator. 

To avoid complications related to the well-known difficulties in treating fermions on the 
lattice, in this study we consider the quenched approximation of QCD, i.e., the pure-gauge 
theory case. We hope to return in a future paper on the inclusion of fermionic effects, which 
may be more important than usually expected for soft high-energy processes (see Refs. |28ll56j l. 

A good choice is to use the anisotropic action discussed previously, Eq. (j2.5p . taking the 
anisotropy parameters to be such that the long sides of the Wilson loops lie in a lattice plane 
at 45° from two of the lattice axes, and are of fixed length. This amounts to set 

A4(g,r)= ^ Ai(g,r)= \ A2(0 ,t) = A3(0,t) = i, (3.5) 

V 2r cos 2 V 2r sm 2 


where T = T/Tq with Tq some fixed length, and 9 is restricted to 0 G (0, tt) without loss of 
generality [TB]. This yields for the plaquette coefficients 


Cai{9,T) r2sin0’ 

C42(^,T) = C43(.,T) = ^^ 


Ci3(0,r) 


(3.6) 


Notice that the following relations hold, 

= A?(0,T) = g^|£|, jr(0,r) = r2sin0 = C23(0,r). (3.7) 

The action dehned by Eq. (123]), with anisotropy parameters Eq. (|3.5I) . will be denoted by 
S'[[/; 9, T], and the corresponding expectation value will be denoted by (.. .)0 f. 

The lattice Wilson loops are defined as 

wl?’ = Atr {IV+iJ+ir-tif-t} , (3.8) 

where the “tilted” Wilson lines are defined as (see Fig. [2|) 

to —1 

Wt= W i7«(iui,2 + 4)> (3-9) 

j=-to 

where ui ^2 = (±1)0,0,1), to = with to G N, and = df /a denotes the transverse position 
in lattice units, see Eq. (13.21) . while are the appropriate Wilson lines made of the usual link 
variables in the transverse plane, closing the loopscj It is clear that \/2to = ^ is the distance in 

^®One can properly choose Ao.s and use “tilted” links also in the transverse plane. This would however leave 
the discussion and the conclusions of this Section unchanged. 
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Xi 


Figure 2: The “tilted” lattice Wilson loops and Eqs. p.Sp and p.9p . projected on 

the longitudinal plane. 


lattice units between the center of a long side of the loop and its endpoints, i.e., loosely speaking, 
the half-length in lattice units of the Wilson loops. The “tilted links” are appropriate 

functionals n] of the lattice links, which in the continuum limit have to satisfjo 

Eq. ()2.4p ] 

U^^\n) = 1 -|- iaTV^ [cos 1 ^ 4 (x(n)) -|- sin |^i(x(n))] + 0{a^) , 

= 1-1- iaTy/2 [cos ^Ai{x{n)) — sin |^i(x(n))] -|- O(a^), 
and which under a gauge transformation behave as 


see 


(3.10) 


G{n)U^^\n)G^{n + l + i), 
(n) G(n)Z^(2) (n)^^(n + 4 - i). 


(3.11) 


The simplest possibility in building is obviously to use a combination of the gauge 

transporters along the two shortest paths connecting opposite corners of an elementary plaquette, 
namely 

U[^\n) = Ui{n)Ui{n + 4), U^^\n) = Ui{n)U4{n + i), 

U[^\n) = U 4 in)ul{n + 4 - i), U^^\n) = C/|(n - i)U 4 {n - i). 

It is convenient to adopt a definition ofZ^^-^^ which is symmetric under the exchange Ui\n) -H- 


(3.12) 


^The factor in front of the square brackets takes into account that the diagonal of a plaquette in the longitudinal 
plane has length yaj+of = \/2aT. Notice that we are using path-ordered Wilson loops, as it is customary on the 
lattice, rather than the time-ordered Wilson loops appearing in the formulae for the scattering amplitudes (see, 
e.g.. Refs. |271153p . This has no consequence on the results, as the theory is invariant under charge conjugation, 
and so under reversing the loop orientation. 
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Fignie 3: The “tilted” links of Eq. (j3.13p . built from the shortest paths connecting opposite 
corners of a plaquette. 



(n). A viable choice is (see Fig. [3]) 


= Proj5^(^^) 


U[^\n) + U^^\n) 


(3.13) 


with Proj 5 [/( 7 v^) denoting the projection on SU{Nc). This symmetry requirement comes out 
naturally if we want that the Wilson loop correlator satisfies on the lattice the same “crossing 
property” [181 [19] that it satisfies in the continuum. It is easy to show that in the continuum 

the correlation function of the two Wilson loops W) 2 i defined in Eq. (|3.2|) , at angle ir — 6 is 

(T) 

equal to the correlation function of Wj 2 angle 9 but with the orientation of one of the loops 
being reversed. In formulae, 


.CTi.. .m. 


(3.14) 


In order to impose this symmetry on the lattice, let us first notice that the anisotropic lattice 
action defined by Eqs. (12.50 and (|3.5I) is invariant under the transformation U = SU^ acting on 
the links, defined by 

C/4(n) = C/f(n“), Ui{n) = Ufin-), C/ 2,3 (n) = C/ 2 “ 3 (n“), 

= = = (3.15) 

114 — ni, n4 — n 4 , n2 3 — 71-2,3 > 

if at the same time one also sends 9 —>■ n — 9. Indeed, it suffices to verify that 6 * 42 (vr — 
9,T) = Ci2{9,T) and C' 4 i( 7 r — 9,T) = C4i{9,T) [see Eq. (|3.6I) ]. Consequently, the one-loop 
corrections JC^u will transform in the same way as as can be also verified explicitly. We 
have then S[SU;9,T] = S[U',tt — 9,T], and since the integration measure is clearly invariant, 
the expectation value of some observable 0[U] satisfies {0[U])^_q f = {0[SU])q f. In order to 
maintain the “crossing property” also on the lattice, the “tilted links” must therefore transform 
as 

C^(i)[HC/;n-] =C/(^)[C/;n] , U^‘^\SU;n~] =U^‘^^^[U-,n - I + 1]. (3.16) 

One can then readily show that the definition Eq. (I3.13[l satisfies the properties Eq. (|3.10|) . 
Eq. (|3.11l) and Eq. (13.1611 . In Appendix [B] we show that using Eq. (|3.13p in the case of the 
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compact U{1) gauge theory one correctly recovers the continuum result of Ref. [T7j in the weak- 
coupling limit. 

One can then define the relevant Euclidean correlator as the continuum limit of the appro¬ 
priate lattice correlator, 


gE{9,T = Tof)= lim gLie,To,f-a,V). 
a—^0,V^oo 


gL{e,To,T;a,V) = 


(W 


Ll ie,T 


^Ll'^)d,T^L2’)e,T 


{To). 


- 1 , 


(3.17) 


where V is the lattice volume, and we have dropped the dependence on the impact parameter 
and on the dipole variables, since they play no role in the following. 


3.3 Analytic continuation 

We now argue that gE{w,T) is analytic in a complex domain V which makes the analytic 
continuation relations Eq. (13.3p meaningful. Here w and T are now complex variables, which we 
parameterise bs w = 9 — ix, with real 9, y, and T = TqT = ro|T|e*^, with ip G (— 27r, 27r]. Since 
one has to take two possibly dangerous limits, i.e., the infinite-volume limit and the continuum 
limit, which currently are not under full theoretical control, onr argument is not rigourous. 
Nevertheless, a few reasonable technical assumptions are snfficient to complete the proof. 

The first thing to check is that the couplings, and the plaquette coefficients, 

are analytic functions of w and T = |T'|e*' 2 . This is obvious at tree level, since 
= 2Nc/g^ and the only singular points of are w = nir with n G Z, and T = 0. 
Analyticity of the one-loop corrections IC^i,{w,T), and so of at the one-loop level, is 

studied in Appendix O 

The next step is to require that the theory has the desired continuum limit. This requires 
positivity of the real part of the action to guarantee convergence. The tree-level convergence 
conditions have been discussed in Ref. [I3], and read 

ReCfj,u{w,T) > 0 (3.18) 

These conditions define a complex domain D which has been fully worked out in Ref. m- 
Although its detailed form will not be used here, it is worth mentioning that D is defined only in 
terms of the complex angle w and of if, i.e., |T| is not restricted (except for asking ITI ^ 0). The 
Euclidean region corresponds to 0 G (0, tt), y = 0, (/? = 0. The Minkowskian region 0 = 0, y > 0, 
ip = TT lies at the boundary of P, and so does also the “crossed” Minkowskian region 9 = n, 
y < 0, (^ = vr; we will refer to these as the “physical” boundaries of P. Notice that both in the 
Euclidean and in the Minkowskian regions the restrictions on the angular variables do not lead 
to any loss of information m- As it is shown in details in AppendixlAj the one-loop corrections 
^ijlv{w,T) are analytic in P. For small enough lattice spacing, the one-loop corrections will 
therefore not spoil the positivity of the real part of the action enforced at tree level, for any 
choice of parameters in a compact subdomain of P. 

At finite volume and finite lattice spacing, and at one-loop accuracy, we have therefore proved 
that the relevant correlators are analytic functions in a domain P, within which positivity of 
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the real part of the action is guaranteed. This domain of analyticity will survive the infinite 
volume limit if the convergence is uniform. Proving this is currently out of reach. However, if a 
lattice system has short-range interactions, then correlation functions of operators localised in 
some finite region TZ of spacetime will become insensitive to the lattice size when this exceeds 
the size of 7^ by a few correlation lengths. Notice that Tq is fixed, so that our operators are 
indeed localised. If interactions remain short-ranged throughout P, then it is enough to take 
the lattice size required by the largest correlation length (within some compact subdomain of 
P) to make finite-size corrections uniformly negligible. This essentially amounts to assuming 
that the theory remains confining as one moves in P. Although we cannot prove this, we find it 
plausible: for example, it is easy to see that it is true at strong coupling by means of a character 
expansion. 

At this point one has to take the continuum limit. This limit is expected to exist and be 
finite within P (again, a rigorous proof is out of question). In particular, Wilson-loop operators 
renormalise multiplicatively [571158], so that the normalised correlation function appearing in 
Eq. (I3.17jl does not require any further renormalisation on top of the renormalisation of the cou¬ 
plings in the action, discussed in the previous Section. A rigorous proof of uniform convergence 
is currently out of reach; however, deviations from the continuum limit are expected to be of 
order 0{a), independently of w and T, and in this case it is possible to make them uniformly 
negligible. 

The conclusion, within the present accuracy, is that is analytic in the complex domain 
P, which, as shown in Ref. m, is sufficiently wide to make the analytic continuation relation 
Eq. (I3.3jl and the crossing-symmetry relations Eq. (|3.4h fully meaningful. 

As it was implicit in the discussion above, singularities in the correlator may develop at 
the boundaries of P. As the analytic continuation Eq. (13.3h is formally equivalent to the usual 
definition of the Minkowskian correlator making use of the ie” prescription |18] . no singular¬ 
ities are expected at the “physical” boundaries of the domain. The anisotropic action itself is 
singular at 0 = 0,7r, x = 0 and 0 = 0, tt, x = cc, but as this is an artifact of the construction 
it is not clear if true singularities are present there. At finite T (i.e., for Wilson loops of finite 
physical length), no singularity is expected in the Euclidean correlator (x = 0) also at 9 = 0,7r; 
however, as T —>• oo, a true singularity is expected to appear there, which has its physical origin 
in the relation between the correlator Eq. (|3.1jl at 0 = 0, tt and the static dipole-dipole poten¬ 
tial [591 Enunn Eg. This is also supported by numerical results [531 ED- On the other hand, 
the points 9 = 0, x = oo and 9 = tt, x = —oo, in the limit T oo, are the ones actually rele¬ 
vant to soft scattering at asymptotically high energy, where the approach initiated by Ref. [20] 
applies. A better understanding of the correlator near these points would help in the study 
of the asymptotic high-energy behaviour of scattering amplitudes and total cross sections. In 
particular, in order to establish that the expressions for the scattering amplitudes derived in this 
approach satisfy unitarity, it is crucial to show that for vanishing 9 and large x the correlator is 
a properly bounded function of the impact parameter and of the dipole variables. Furthermore, 
the existence (or not) of the strict x ^ oo limit at fixed impact parameter, and the properties 
of the correlator in this limit, are closely connected to the issue of universality of hadronic total 
cross sections observed in experiments (see, e.g.. Refs. [63( 164] and references therein). For more 
details on these problems, we invite the interested reader to confer Ref. |28j . 

Other singularities could appear when |T| —)• 0 or |T| —oo. Working at fixed Tq, this 
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corresponds to |T| —)• 0 or |T| —)• oo, which are again singular points of the anisotropic action. 
However, since the analytic continuation to the “physical” boundaries requires only the phase 
of T to be changed, one can take as well Tq = |T| and T = exp{i^}, and study the two 
limits above by changing the length of the “tilted” Wilson loops@ The above limits therefore 
correspond to the limit of “tilted” Wilson loops of vanishing or infinite length. In the first case 
no singularity is expected; in any case this limit is irrelevant for our purposes. On the other 
hand, the limit of infinite length is the one entering the physical scattering amplitudes. In this 
case, the short-ranged nature of strong interactions (which is assumed to remain unchanged 
throughout the analyticity domain) implies that distant parts of the two Wilson loops do not 
“feel” each other, i.e., those parts of the loops that he beyond a certain distance from the centers 
interact mutually only very weakly, and essentially contribute only to the self-interaction of the 
loops. These contributions are cancelled by the normalisation factors, so that the correlator 
becomes basically insensitive to the length of the loops beyond some critical value, and a finite 
limit |r| —oo is therefore expected. In the Euclidean case, this has already been checked on the 
lattice, although in an isotropic setting [53]. As discussed in Ref. m, the boundedness and the 
analyticity properties of the correlator as a function of T imply through the Phragmen-Lindelof 
theorem (see, e.g.. Ref. [65] 1 that the analytic continuation to Minkowski spacetime and the 
infinite-length limit commute. Setting Ce,m = liinr^oo this means that one can obtain 

the physical correlator by means of an analytic continuation in the angular variable only, i.e., 

Cm(x; z±-,Ri±, fi',R2±,f2) = CEi-ix', z±;Ri±,fi;R2±,f2) ■ (3.19) 

The analyticity domain for (tc = 0 — ix), already discussed in Ref. [I3| , is clearly not changed 
by one-loop corrections, and it is simply the strip 0 € (0, vr), y G R, shown in Fig. Sj 


4 Longitudinally rescaled action 

The results of Section [2] can be used to obtain some insight in the approach to high-energy 
scattering based on longitudinally rescaled actions [6l|7l[8l[9l[ini[IIl[l2]. The physical idea behind 
this approach is that in high-energy scattering processes the longitudinal directions appear highly 
Lorentz-contracted, so that it should be possible to achieve an effective description through an 
appropriately rescaled action. While initially only a classical rescaling was considered PEIEIIS], 
in recent years the effects of quantum corrections have been computed by means of anisotropic 
renormalisation in the continuum theory unittiiiii]. Here we will consider the same problem 
in the lattice approach, which will allow us to clarify, to some extent, the structure of the action 
in the limit of large anisotropy. Notice that the anisotropy class (2-1-2) is the same considered 
in Ref. |5]. 

On the lattice, the tree-level anisotropic action is given by Eq. m with the following 
anisotropy parameters, 

(0 = Af (e) = e, (0 = A^^) (e) = 1. (4.1) 

^®In the continuum limit the choice of To should be irrelevant, as long as it is compensated by the appropriate 
choice of T. 
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Figure 4: Analyticity domain of the Wilson-loop correlator with the infrared cutoff removed, 

[see Eq. (|3.19l) ]. The solid black lines indicate the boundaries of the domain, and crosses signal 
the singularities. 


In the following, the superscript LR is used to specify that this particular choice has been made. 
We will refer to directions 4 and 1 as longitudinal, and directions 2 and 3 as transverse, and 
use the notation ny = (n 4 ,ni), n± = (n 2 ,n 3 ), ay = 04 = oi = a/^, a± = 02 = 03 = a. The 
plaquette coefficients in the anisotropic action read 

^(0 = ^ , c^^4(0 = (0 = (0 = (0 = 1 , (0 = . (4.2) 

The interesting case is that of large 0 Taking naively the limit ^ 00 in the tree-level action, the 

transverse-transverse plaquette term drops from the action, while the longitudinal-longitudinal 
term yields essentially a “delta function” forcing the longitudinal links to be trivial. The resulting 
effective action would read 

g,,„ + 

n_L 

n|| q=4,1 

which describes a set of independent 2D principal chiral models involving the transverse link 
variables, each one living in the longitudinal plane at a given point n± in the transverse plane. 
Here A+ has been redefined by omitting the Aq, factor [see Eq. (|2.15p ]. Taking into account 
quantum corrections, however, a different coupling has to be used for each of the three different 
kinds of plaquette terms, namely for the longitudinal-longitudinal term, = 

1 ^ 23 ^^ for the transverse-transverse term and = ^[ 2 ^^ = for the 
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longitudinal-transverse terms. Recall that the quantum corrections are of the form 

64' 




3 (47r) 


-7 + 
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fii' 


(4.4) 


= -/3o log c2 + (0 + (0 , 


liu 


with and AZ^^^ containing respectively the contributions of the G- and 21-integrals, 

Eqs. (|A.ljl and (jA.3jl . Obviously, AG^^^'^ = AG^^^^ = = A^jg^^ and similarly for 


AZ\^^\ Using the large-^ behaviour of these integrals, derived in Appendix O one gets 


o{LR) 

P L L 


(o,0 


2N, 

2N, 

2Nr 


= Po log 
= /3o log 
= Po log 


1 


(aAc)2 

1 

(aAc)2 

1 




±± 


“(0 + AZ 


±± 






(4.5) 


A. 1 




{LR),fin 


(0 + AZ| 


{LR)^fin 

lU 


(0, 


where the superscript fin on a quantity indicates that it is finite in the limit ^ oo. If we 
keep the transverse spacing a_L = a fixed, then taking ^ oo means taking ay = to zero, 
i.e., taking the continuum limit in the longitudinal plane only. One then sees that in general 
it is not allowed to discard the transverse-transverse plaquette term, since = 

Sni| contains the right power of ay to become the two-dimensional integral over 

the longitudinal plane in the limit ay —)• 0. 

The action can now be recast in a form appropriate for a set of coupled two-dimensional 
principal chiral models. To this end, it is convenient to introduce the following couplings, 




l3{2D){a\\,a±) — 


= 1 
0(^0) 2 47r ayA(2^)(a_L) ’ 


/3(2D)(a|p «±) — 


/3(2Z))(o±) — 


2N, 


N? -11 


0(50) 


■log 


1 


2Nc vr ay A(2-D)(a_L) 


(4.6) 


2Nr_ 


0(50) 


= Po log 


(a_LAc) 


2 ’ 


where the a_L-dependent scales and are given in terms of the original A-scale as 

follows, 

A(2^)(a^) = , 

A(2^)(a_L) = Ac(a_LAc)TU e ^ . 

The action can be equivalently written as follows, 

5iat = 4'^(nA) + S^^\n2_) + 5inti(nx) + , (4.8) 


nj_ 
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where 


correspond to principal chiral models, 


= /3^2D){a\\,a±)Y^ ^ tr {[A+C/^(n)][A+C/^(n)]t} , (4.9) 

ny 0=4,1 


and the interaction terms read 

5’inti(?^±) = P{2D){a\\,a±)2^—V23{n) , 

nil 

*S'int2(?^±) = /5(2 D)(o||) 0±) X] Z] Z] [2A^c7^;.o(n)-tr{[A+[/^(n)][A+[/^(n)]t} 

/i=2,3 ^T-y 01=4,1 

+ P{2D){a±)'^‘^Nc-^V4i{n). 


(4.10) 


The only approximation made here is to discard o(i^*^) terms in the couplings, so that this is just 
a rewriting of the original action in the limit of large Nevertheless, this expression displays a 
remarkable feature: the coupling /3(2D) is precisely the one appropriate for a 2D principal chiral 
model with lattice spacing ay, to one-loop accuracy (see, e.g., Ref. [66]). The principal chiral 
models are clearly not independent, with the precise form of the interaction dictated by the full 
4D action. Notice that identifying the longitudinal links with = exp{iayg^} and expanding 
in powers of ay, the summands in the interaction term S'int 2 are of order 0(a|), as appropriate 
to obtain an integral over the longitudinal plane in the naive ay —?■ 0 limit, so there is no reason 
to discard these contributions El It is not surprising that the interaction terms cannot be be 
neglected a priori: after all, no matter how anisotropic the lattice is made, by construction 
the action has to describe QCD in the continuum limit. The possibility or not to neglect the 
interaction terms will depend on the properties of the specihc observables relevant to the study 
of high-energy processes. 

The expectation values of the different plaquette terms can be used to estimate the range of 
applicability of the expressions above. Using Eq. (I2.56P one gets to lowest order [see Eq. ()A.37p ] 


(^4l) 

(^42) 

{V 23 ) 


—2-y-~ ^ - 


= g 


= g 


Nc 

2 A^c -1 

2Nc 

t2 


r{LR) 


(0 + 


4"«)(0' 


- g 


2^^c -1 

2Nr. 


-ZlO 




= g 


Nr. 


Nr 47r 


(4.11) 


with ziQ a constant dehned in Eq. (IA.21D . so that in order to have small fluctuations one needs 
g^log^ <C 1. Together with the basic assumption 1, and the fact that we work here at 

notice that working in the axial gauge Ui = 1 and expanding 5int2 to D(ay), the resulting expression is 
quadratic in 54 and the corresponding integration can be carried out. This leads to the appearence of complicated, 
non-local interaction terms involving the transverse link variables. 
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^ S> 1, the requirement (P 23 ) ^ 1 defines the range of applicability of perturbation theory, 
which in terms of the lattice spacings reads 

^ AtvPqNc 

1 >o_LA>a||A> (auA) > {a±Kf. (4.12) 

The important fact is that Eq. (14.121) does not allow us to strictly take the continuum limit 
in the longitudinal plane before taking o_l to zero. This was already suggested in Ref. m, 
although there it is claimed that perturbation theory makes sense only for ^ slightly larger than 
1; according to our results, a much larger region seems to be accessible. 

A comparison of our results with those of Refs. [laiiiKii] is not straightforward. First 
of all, since we use a different regularisation, we expect different finite contributions to the 
renormalisation of the couplings in the limit a —>■ 0 (at fixed ^); ultraviolet divergences, on 
the other hand, have to be the same. Indeed, to account for a change in the cutoff, Orland 
and collaborators integrate over an anisotropic ellipsoidal shell in momentum space, while on 
the lattice a change in the cutoff requires us to integrate over an anisotropic parallelepipedal 
shell. It would be interesting to compare the divergent terms in the limit 00 , but in 

Refs. [TOl HTJ [12] only the case ^ > 1 is studied. 

We conclude by noticing that a similar recasting of the action can be done also in the case 
discussed in Section [3j considering the limit of large T. The results are briefly discussed in 
Appendix O 


5 Conclusions 

In this paper we have performed the renormalisation of SU{Nc) gauge theories on a general 
four-dimensional anisotropic lattice, with different lattice spacings in the four directions, using 
perturbation theory to one-loop order and the background field method on the lattice (Section 
n. To avoid the complications related to the introduction of fermions on the lattice, we have 
discussed here the pure-gauge case only. For general anisotropy, the various couplings in the 
gauge action need to be properly tuned in order to recover 0(4) invariance in the continuum 
limit, as already observed in Ref. [5|. In practice, however, only two parameters need to be 
tuned for this purpose, which reduce to one if there is at least a pair of equal lattice spacings 
(and to none in the 3-|-l case). A simple nonperturbative scheme for this tuning, based on the 
string tensions obtained in different lattice planes, has also been proposed. 

In Section |3l the possibility to vary continuously the anisotropy parameters has been ex¬ 
ploited in the context of the nonperturbative approach to soft high-energy hadron-hadron scat¬ 
tering based on Wilson loops [SgEIlEaESllMllSSlEelETl, in order to refine the arguments 
of Ref. |13j on the analyticity properties of the relevant Wilson-loop correlators. The results 
reported here give further support to the possibility of performing the desired analytic con¬ 
tinuation between Euclidean and Minkowski space, and thus on the very possibility of using 
Euclidean techniques to study soft high-energy processes. This is particularly important in the 
light of recent progress on the problem of hadronic total cross sections [281 [56] , which is based 
on the possibility of recovering the physical amplitudes starting from Euclidean space. 
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In Section [J] we have applied our results to the longitudinally rescaled actions considered 
in Refs. p El El nnni EH 1^ to study high-energy scattering in QCD. At the classical level, 
in the limit of large anisotropy the action reduces to that of a set of coupled two-dimensional 
principal chiral models, living in the longitudinal plane at each point of the transverse plane. 
Our main result in this context is that this interpretation holds also at the one-loop level, as 
the bare coupling resulting in the free part of each principal chiral model behaves appropriately 
as a function of the longitudinal lattice spacing. The precise form of the interactions among 
the principal chiral models is dictated by the full gauge action. However, the limit of large 
anisotropy cannot be taken independently of the continuum limit, at least in the perturbative 
approach. Indeed, the requirement of small gauge field fluctuations defines a range of validity 
of the form 1 ;§> a^A ^ ayA ;:|> (auA)^"'''’' for the longitudinal and transverse lattice spacings 
ay and a_L, where A is the QCD scale and 7 > 0. Nevertheless, our findings suggest that there 
may be a deeper relation between gauge theories and principal chiral models than just at the 
classical level. 

There are several open directions for future studies. An obvious possibility is the inclusion 
of fermions in the analysis. This is particularly relevant to the nonperturbative approach to soft 
high-energy scattering, since the presence or not of dynamical fermions seems to have large effects 
on total cross sections [28l [56]. It would be interesting to extend the perturbative analysis to 
non-orthogonal lattices, which would allow us to use on-axis Wilson loops in the relevant lattice 
correlator. However, in this case more terms appear in the action, so leading to a more intricate 
calculation. 


Acknowledgments 

I want to thank E. Meggiolaro and T. G. Kovacs for many useful remarks. I also acknowledge 
interesting discussions and correspondence with G. Burgio, E. Eollana, S. Katz, and P. Orland. 
This work is supported by the Hungarian Academy of Sciences under “Lendiilet” grant No. 
LP2011-011. 


26 


A The Q- and ^-integrals 


In the expression for the one-loop contributions Eq. (j2.53j) . there appear a few integrals 
involving the modified Bessel functions of the first kind Iniz), which are special cases of the 
integrals 


Gn{X) = 


_OL=l 


f 


Zni\) = dp 


dpp llKlt‘^\2Xlp) 
=1 
4 

Aq,Io(2A^p) 

a=l 

n xjt‘^\2\ip) 


dpp 

4 


n / ( 0 , 0 , 0 , 0 ) 

1 1 


-e(p-i) 


(dvrp)^ 


n= (0,0,0,0), (A.l) 


.Ci=l 


defined for a general four-vector of integers n, where 

In{z) = e-^In{z) , lt\z) = {-d/dzrin{z) , 


(A. 2 ) 


and ©(z) is the step function. In particular, in Eq. (I2.53P we have denoted as follows the relevant 
cases, 


Q^iv — ^nL I i 

*lt(yi - ^OLpL~^OLl/ 


Gu = Gn 


G = G 


n\na=0 ) 


, 


Z — Zn\na=0 • 

These integrals are not all independent; in particular, the following sum rules hold, 
4 4 4 


(A.3) 


/i=l fi=l fi=l 

2Xlg, + Y,2Xlg^, = Y,‘^xlg,. 


E ^ ’ 


(A.4) 






A global rescaling Aq —)• CAa (C > 0) of the anisotropy parameters can be essentially reabsorbed 
in Zn and in by changing variables to p' = which brings about a multiplicative factor 
for Zn, and an additive contribution proportional to logC to ^o- More precisely, we find 


Zn{CX) = (^Zn{X), en(CA) = 


0 n(A), n/ 0 , 

J— lncr ^2 

(dvr 


^o(A)-b 773^1ogC^ n = 0. 


(A.5) 


A.l Analyticity properties 

We discuss now the analyticity properties of K,^y. It is clear that for A^ 7 b oVa these depend 
only on the analyticity properties of the g- and 2i-integrals defined in Eq. ()A.4p . Since these are 
integrals of analytic functions of p and {Aq,}, it suffices to show that they converge uniformly 
in {Aq,} within some complex domain. In turn, a sufficient condition for this is that we can 
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bound the modulus of the integrand uniformly in {Aq} by some function /, whose integral is 
also convergent. To do this, we need the following inequalities, 

\Io{z)\<Io{Rez), \Ii{z)\<Io{Rez), |7o(^)| < 1 ifRez>0, (A. 6 ) 

which are easily proved using the integral representation for /„(z). We also need the monotonic¬ 
ity property 

^7o(x)<0, VxGM, (A.7) 

and the asymptotic behaviour of Io{z), 

(a.8) 

valid for | argzj < vr (see, e.g.. Ref. [67!). 

1 The quantities Z and Z^ are given by the product of the analytic factor Aq, and 

an integral of the product of functions Iq and, in the case of Z^, also Iq — Ii, so that we may 
write 

^({A.}) = . (A.9) 

For {Aa} such that for every a one has Re A^ € [uq,, Va], with Uq, G M, 0 < Mq, < < oo, the 

first two inequalities in Eq. ()A. 6 I) and the monotonicity property Eq. (|A.7h tell us that a possible 
choice for f{p) to bound the modulus of the integrands both in Z and Z^ is f{p) = 2 faip), 
fa{p) = Io{‘^'UaP)- In particular, this shows that Z and Z^ are analytic functions of {A^}. 

2 To study G we split the integral into two parts, Eor the first piece, the third 

inequality in Eq. (lA. 6 h indicates that we can take f(p) = p. The integrand of the second piece 
is conveniently written as 

P Aq,7o(2Ao,p) - ^ ’ (A.IO) 

where / is analytic VA and p / 0 , and furthermore it is certainly bounded for ReA^ G [uaiVa] 
and p G [0, oo), since it has a finite limit as p —>■ oo, see Eq. (jA.Sp . In this case we can then take 
f{p) = M/( 47 rp)^ for a properly chosen constant M. 


3 Einally, analyticity properties of and Gfn^ are inherited from Z and Indeed, since 
one can bring derivatives under the sign of integral due to uniform convergence, one shows 
immediately that 


x>.^^z{x) = z{\)-4xlg^{x), 
X,^^Z^{X) = Z^{X)-AXlG^,iX) {v^p). 


(A.ll) 
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Notice that and Q^u are of the form 


= J-^G^..{{Xl }), (A.12) 

with and G^iv analytic in {A^}. 

In conclusion, are analytic in any compact domain with Re > 0, Va. For our purposes, 
it is convenient to extend further the domain of analyticity. To this end, notice that for real 
positive Aq,, one can rewrite Z, and G^xv as follows by exploiting their behaviour under 

global rescaling, Eq. (IA.5I) . 

Z({A„}) = Z{{J\l }), Z.aXx,}) = MiJXl }), ^ 

G,{{X^}) = JG,{{JXl}), G,A{K}) = JG^u{{JXI}), 

where Eqs. ()A.9p and (IA.12(1 have been used. The domain of analyticity of these quantities can 
thus be straightforwardly extended to Re {JX^) > 0. Eurthermore, for real positive Aq, one has 


G 




0(p-l) 


1 


1 

(47r)2 


log 


(A.14) 


where we have used Eq. (IA.5I1 again. By the same token used above in point 2, the first term in 
Eq. (]A.14I1 is analytic for Re^A^ > 0. The logarithmic term is an analytic function in the cut 
complex plane for | argj7| < tt , so we conclude that fC^^u are analytic also in the domain defined 
by Re {JX\) > 0, | arg J\ < tt. 

We now analyse the specific case discussed in Section [3l corresponding to the following choice 
of anisotropy parameters, 

A4(g,f)= ^ Ai(0,r) = -J-^, A2(0,T) = A3(0,r) = l, (A.15) 

V 2T cos I v2Tsm| 


which, in the light of the extension of the analyticity domain discussed above, can be recast 
more conveniently as follows, 


j{e,f)xl{e,f) = C 42 ( 0 ,r), j{e,f)xl{e,f) = Ci2(0,f), 

J{e,f) = J{9,f)XU9,f) = J{9,T)XU9,f) = C23{9,T) . 


As functions of complex angle and length, Cfj,i,{w,T) are analytic everywhere, except at w = 
nvr with n € Z, and |T| = 0. Since the domain D considered in Section [3] is defined by 
^eCfxu{w,T) > 0, in P one has that \ avg J'{w,T)\ < ^ and Ke{J'{w,T)Xa{w,T)) > 0, so that 
the G- and Z-integrals are analytic there, and in conclusion the one-loop corrections ICfxu{w,T) 
are analytic in P. 


A.2 Large-T behaviour 

We now determine, for real T, the large-T behaviour of the Z- and ^-integrals for the choice 
of anisotropy parameters of Eqs. (13.5p and (jA.15h . To this end, it is convenient to define the 
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following auxiliary quantities, 


roo 

Bn{e,f)= / dp 

Jo 

roo 

{e,f) = / dpp 

Jo 


D. 


\{KInA‘^\lp) 

a 

tf AQ,7n„(2A^/)) - 0(p- 1)7T— 


(47rp)2 

where {Aq} are chosen according to Eq. (|3.5I) . It is straightforward to show that 


^ JJn\nc,=0 1 Bn\na=0 JJri\na=Saij. i 

G ~ -^n|nQ .=0 ) Gfi ~ -^n|nc ,=0 ^ri\na=&aii i 


Gjj.u — JJn\nc,=5c,ij,+5av 4 " JJri\na=0 JJri\nc,=&c,ij, ^n\nc,=5c,v 


A rather simple calculation shows that at large T 

1 f 1 
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sin 6 1 dvr 
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»1 


bn{6) = InA^)Ini{Q) / dp In 2 {‘ip)Ini{‘ip) 

Jo 


4 ‘ 
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Dn{e,f)= -^logr2+d„(0) + o(fO) 
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dp 
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sinfy Vcos^i; ^sin^l 
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dvrp 


(A.17) 


(A.18) 


(A.19) 


It is now straightforward to obtain the large-T behaviour of the relevant quantities. For the 
2^-integrals we have 


z{e,T) = 
z^{e,f) = 

Zi{e,f) = 


1 1 


sin 6 1 dTT 


log + Zoo + zqo{6) + o(T^) 


log + Zoo + ^ 10 (6') + o(T°) 


sin 0 1 dvr 

1 f 1 


sin 9 1 dTT 


log + Zoo + ^01 (0) + o(r°) 


Z 2 ( 0 ,r) = Z3{9,f) = 
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r2 sin e 


Zio + o(T°) 


(A.20) 
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where we have introduced the following quantities, 


^nm — 


( 2 p)/<”> ( 2 rt , 

Jo 
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,2 0 0 
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,2 e 


- p) 


(A.21) 


For the ^-integrals we find 


0(6,T) = - logT^ + -gooiO) + o(f0), g,i{e,f) = gu{9) + oif^) . 

(47r)^ 


e4(0,r) = 5ioW + o(fO), 
where we have introduced the following quantities, 


gi(0,f) = 5oi(0) + o(T°) 


(A.22) 
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(A.23) 


sin 


2 0 


-0 (p-l) 


dvrp r ’ 


while the remaining integrals are all o{T^). One can now easily determine the contributions of 
both kinds of terms to JC^^, namely 


Ag^,{e,f) = N, 

= /3olog^ + Ag/-( 0 ,r), 
where Pq is defined in Eq. (j2.52l) . and 

Ag(l^{e,f) = N, 

Ag(l^{e,f) = ( 0 ,f) = N, 

Ag(i^{9,f) = Ag(i^{9,f) = N, 
Agts^i9,T) = o{f^), 


lg,.{9,f) - ^{g^{e,f) + g,{9,T)) + ^g{9,f) 


(A.2d) 


Y5oo(6') - ^{m{9) + goi{9)) + pii(6») 


+ o{f^), 


^5oo(0) - %{9) 


3^ 

^500(0) - poi( 0 ) 


+ o(r 0 ), 
+ o(r 0 ), 


(A.25) 
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where can be split into a divergent and a finite part, 

AZ4i(0, f) = f) + AZ{i^(0, f ), 

AZ42(0, f) = AZ43(0, f) = AZflie, f) + AZ{f (0, f ), 
AZi2(0, T) = AZi 3(0, T) = AZfr(0, T) + AZ{f (0, T), 
AZ23{e,f) = Azi^{e,f) = o(f0), 
with divergent parts given by 


(A.27) 


— 1 9 1 
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IjSc Sin u 47r 

Azfr( 0 ,r) = cot^^^iogr 2 , 


AZfr(0,r)=tan^^^logr2, 


and finite parts given by 
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+ o(r°), 
+ o{f^) , 


(A.28) 


(A.29) 


from which one can easily reconstrnct the behaviour of lC^y{9^T) up to o(T^). 

The results above allow us to easily derive the large-^ behaviour of the couplings when the 
anisotropy parameters are chosen appropriately for the longitudinally rescaled action of Section 
m i.e., ^ and = Ag^'^^ = 1, see Eq. (14.11) . This is accomplished through 

the following steps. First of all, notice that A^ are just a particular case of X^{9,T)^ namely 
’ |)- Next, it is straightforward to show that 


B, 
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(47r)2 


log^2 + T>n(f,6> 


where = np_ with {/!} = {1,2,3,4} = (3,4,1,2}. Finally, one easily shows that 

_ 1 

[A;.(f,|)F [AA(f,0P' 

Putting these results together one finds that 

AZ<“>({) = A2„„(§,i) = AZ^c(f,0, 

A<;l“>(a = A5„„(f, f) + Aj, log?" = A<;?“(f, f). 


(A.30) 
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Explicitly, 
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AeS^^(0 = Agjj'^^(0 = o(e°) 

for the contributions AG ^^^^, and 
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for the contributions Az\t/^\ with divergent parts 
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1^ 4 
We also report the values of the 2^-integrals, 

= zi^^\o == ^10+o(^°), 

z[^^\0 = ^ log+ ^00 + 5lo(f) + o(C°) . 


(A.33) 


(A.34) 


(A.35) 


(A.36) 


(A.37) 


B Abelian case 

In this Appendix we compute the Wilson-loop correlator considered in Section [3] in the compact 
[7(1) lattice theory and in the weak-coupling limit. The starting point is the 4D anisotropic 
lattice formulation for the U[l) gauge group, 

^ E - ReU,,{n)) . (B.l) 

n,fi<u 
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Here the plaquettes are built with the ?7(1) links C/^(n) = exp{i(j)fj_{n)}, and can be written 
as ReU^uin) = cos with = (p/iin) + (t)u{n + ft) - + D) - (p^in). The Haar 

measure is simply f dU^{n) = Setting 


U^^\n) = Ui{n)Ui{n + 4) = S^Rn)+4>Rn+l)] ^ ^ 




U^^\n) = Ui{n)Ui{n + i) = e*['^4(n+i)+<Ai(n)] ^ givp«(n) ^ 

Uf\n) = Ui{n)Ul{n + 4 - i) = e*[<^4(’^)-<Ai(n+4-i)] ^ ^ 


Uf\n) = Ul{n - l)Ui{n - i) = e*['^4(n-i)-<Ai(n-i)] ^ (n) ^ 

the analogue of Eq. (|3.13p is here 

$(i)(n)' 


U^^\n) = ( . - . 0) = exp + ip^i\n)) sign cos ■ 

h’T.) _g e*‘^2 ("■)! 


^0 )! 


(B.2) 


(B.3) 


where = $41 (n) and 4>(^)(n) = 4>4i(n — i). The Wilson loops are written as = 

e®^*(TfcTfc, /c = 1,2, with 


— 


1 

2 X] + 4+) + ^'i\jvk + 4+ 

j=-to 


)-<fi {JVk 






(B.4) 


CTfc the product of the sign factors appearing in Eq. (IB.3p . and Tk the contribution from the 
transverse links. 

The calculation is greatly simplified if we take the limit T —>■ 00 firstly Discarding the 
longitudinal-longitudinal plaquette term, enforcing the triviality of the transverse links, and 
using 1 — ReD^Q, = for trivial Ua links, one ends np with 


qU(l) 

“^lat 


T—>-oo 


1 


/i=4,l 


(m) 


Y1 Y1 , 

n|| ,nx a=2,3 


c(") 



c(i) 



(B.5) 


where ny = (n4,ni) and n± = (742,743). The Wilson loops simplify to Wla: e*^'=(Tfc. Since 

there is no interaction between link variables living at different sites of the longitudinal plane, 
and between and Ui variables, one easily sees that the “tilted links” of Eq. ()B.3I) interact with 
each other only if they are separated by at most one lattice spacing in the longitudinal plane, 
which leads to factorisation of the Wilson-loop correlation function and expectation values. 

It is convenient now to rescale the phases as 4(^) ~ with x = en± (notice that 

X is dimensionless), in order to take the weak-coupling limit. One then obtains for the action 


qUil) 

•^lat 


—y 

T—>-oo, e—>-0 


E 


c(m) 


^ ^ ^[da(p^infx)f , 

?T^|| 0=2,3 


(B.6) 


Since there is actually no continuum limit to be taken, in this case the complications of the non-Abelian case 
are absent. 
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and the integration measure in the weak-coupling limit becomes 



(B.7) 


where we have omitted a factor e/(27r) since it gets cancelled in expectation values. We passed 
to the continuum notation for simplicity: as the action is quadratic, the resulting continuum 
Gaussian integrals are fully under control. The propagator is readily obtained, 

Df,u{n \\, my ;x,y) = ; x)^^(m||; y)) = “ v) > (^-8) 

where D{x) is the 2D scalar propagator, 

D{x) = -^log\x\. (B.9) 


From here on angular brackets without subscripts denote the expectation value with respect 
to the action Eq. (IB.61) . In the weak-coupling limit, cosd'^jy = 1 + 0{e'^) and we can neglect 
the sign factors in the expression for the Wilson loops, i.e., WLk Since the action is 

quadratic, one has for the relevant correlation function as e ^ 0 


lim 




T^oo (y^l) 0^T 


g-i{(Ql+02)2> 

g—5{r2j}g—iiQj} 




(B.IO) 


Using now the explicit expression for see Eqs. (IB.2P and (IB.41) . and exploiting the fact 
that the propagator is diagonal in the link directions and in the longitudinal coordinates, a 
straightforward calculation gives 


(D1D2) 


6 

— cot 6 log 

ZTT 


+ % 

y Ri± R2± 

'^-^2 2 

J_ -^2_L 

"T 2 2 

Rl± 1 R2± 
^_L 2 2 


(B.ll) 


which agrees with the known result for in the 4D 17(1) pure gauge theory in the continuum 
limit m- Here we have set /i = /2 = | for convenience, without any loss of information |54] . 


C Large-T limit of the (^, r)-dependent action 

For completeness, in this Appendix we report on the large-T limit of the anisotropic action with 
anisotropy parameters Eq. (13.51) , discussed in Section [3l The idea is that there could be some 
useful simplification if one takes T ^ 00 , corresponding to the limit of loops of infinite length, 
before taking the continuum limit. In full analogy with the discussion of Section 01 in this limit 
the action can be recast as that of a set of interacting principal chiral models, which however live 
now in the transverse plane at every site of the longitudinal plane. This is natural since the limit 
T —>■ 00 corresponds to taking the continuum limit in the transverse plane at fixed longitudinal 
spacing aii = Ta, i.e., the same situation of Section 0] but reversing the roles of the longitudinal 
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and the transverse planes. In the large-T limit, the longitudinal-transverse couplings can be 
rewritten as follows, 


/3(^^(a^,a||,0) = ^C42 = 


/342 


1^43, 


lNr_ 


1 


2Nr 


2Nc 2 dvr a±A^^2(«|| > 


A^^(a||,0) = Ace (^a\\Ac) 


Nc 


/3S(a±,a||,0) = ^Ci2 = 


Pl2 


f^l3 


lNr_ 


1 


2Nr_ 


2Nf^^ 2 47r^°^«^AS(a||,0) 


cot 1-1 


(C.l) 


AW(«|h^) = Ace-^“*2(^eiTM+^^iTM)(a||Ac)^ 

which is precisely the form of the bare coupling as a function of the lattice cutoff a_L = a in 
the two-dimensional SU{Nc) principal chiral model, to one-loop accuracy (see, e.g.. Ref. [66]). 
Here we have neglected o{T^) terms. The remaining couplings read, in the same limit and in 
the same approximation. 


2Nr. 


C 41 = 


1 -I 1 


log- 


1 


1^23 


T‘^sm.9 2Nc vrsin^ a±h.2D{a\\,9) 

27rNc 

N'i 

= f^hD{a\u9). 


= ^hD{a±,a\u9) , 


A 2 d(o||, 6 ') = Ac(a||Ac) ^"-1 e ^"-1 ^ ^ ^ 41 ^ u 


C 23 = T-^ sin 0/30 log ■ 


1 


2Ar,--^ °(a||Ac)2 

The action can be recast as follows, 

S{2D) ^ ^ ^ 5intl(n||) + Sint2(n||) , 

nil 

where correspond to principal chiral models, 

S’5l^)(n||) =/3^^^(a||,a_L,0)^ ^ tr{[A+[/^(n)][A+C/^(n)]1'} , 

nx q:=2,3 


(C.2) 


(C.3) 


(C.4) 


and the mutual interactions are given by the remaining terms, 


5'inti(^ll) = ho (a||,a±,0) y~] ^7^41 (n), 

n_L II 

5'int2(n||)= ^ /3j'^^(a||,a_L,0)^ ^ j2AcP;.a(n) - tr {[A+I7^(n)][A+C//,(n)]t} 
/^=4,l nx 0=2,3 

of 

+ /32D(a±,0)^2A,4iP23(n). 

d I 

n I -L 


(C.5) 


The two-dimensional scales A^^^^ and A 2 D have prescribed values that depend on A, which is 
set in the 4D theory, and on ay, which has to be taken to zero at the end of the calculation. 
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However, the average plaquette terms to lowest order and for large T read [see Eqs. (12.561) and 

(fOOi) ] 


~ g 


logr^ 

Nc dvr sin 0 ’ 


()P23) - 9 


2 - 1 ^10 
Nc f2sin0’ 


, oN^-i e ^ 2^c-i ^ 

()P4±) ^ 9 cot -Zio , {Vi±) ~ 


2Nr 


2Nr 


tan -zio , 


(C. 6 ) 


so that the range of applicability of perturbation theory is limited by g'^ log T <C 1 ; more precisely, 
besides ay 3> a_L one needs ay A <C (a_LA)^“'>' for some 0-dependent 7 , which prevents from taking 
the continuum limit in the transverse plane independently from the longitudinal plane. 
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